Additive and non-additive genetic variance in juvenile Sitka spruce (Picea sitchensis Bong. Carr).

Many quantitative genetic models assume that all genetic variation is additive because of a lack of data with sufficient structure and quality to determine the relative contribution of additive and non-additive variation. Here the fractions of additive (fa) and non-additive (fd) genetic variation were estimated in Sitka spruce for height, bud burst and pilodyn penetration depth. Approximately 1500 offspring were produced in each of three sib families and clonally replicated across three geographically diverse sites. Genotypes from 1525 offspring from all three families were obtained by RADseq, followed by imputation using 1630 loci segregating in all families and mapped using the newly developed linkage map of Sitka spruce. The analyses employed a new approach for estimating fa and fd, which combined all available genotypic and phenotypic data with spatial modelling for each trait and site. The consensus estimate for fa increased with age for height from 0.58 at 2 years to 0.75 at 11 years, with only small overlap in 95% support intervals (I95). The estimated fa for bud burst was 0.83 (I95=[0.78, 0.90]) and 0.84 (I95=[0.77, 0.92]) for pilodyn depth. Overall, there was no evidence of family heterogeneity for height or bud burst, or site heterogeneity for pilodyn depth, and no evidence of inbreeding depression associated with genomic homozygosity, expected if dominance variance was the major component of non-additive variance. The results offer no support for the development of sublines for crossing within the species. The models give new opportunities to assess more accurately the scale of non-additive variation. Supplementary Information: The online version contains supplementary material available at 10.1007/s11295-023-01627-5.

Enhancing genetic disease control by selecting for lower host infectivity and susceptibility.

Infectious diseases have a huge impact on animal health, production and welfare, and human health. Understanding the role of host genetics in disease spread is important for developing disease control strategies that efficiently reduce infection incidence and risk of epidemics. While heritable variation in disease susceptibility has been targeted in livestock breeding, emerging evidence suggests that there is additional genetic variation in host infectivity, but the potential benefits of including infectivity into selection schemes are currently unknown. A Susceptible-Infected-Recovered epidemiological model incorporating polygenic genetic variation in both susceptibility and infectivity was combined with quantitative genetics selection theory to assess the non-linear impact of genetic selection on field measures of epidemic risk and severity. Response to 20 generations of selection was calculated in large simulated populations, exploring schemes differing in accuracy and intensity. Assuming moderate genetic variation in both traits, 50% selection on susceptibility required seven generations to reduce the basic reproductive number R0 from 7.64 to the critical threshold of <1, below which epidemics die out. Adding infectivity in the selection objective accelerated the decline towards R0 < 1, to 3 generations. Our results show that although genetic selection on susceptibility reduces disease risk and prevalence, the additional gain from selection on infectivity accelerates disease eradication and reduces more efficiently the risk of new outbreaks, while it alleviates delays generated by unfavourable correlations. In conclusion, host infectivity was found to be an important trait to target in future genetic studies and breeding schemes, to help reducing the occurrence and impact of epidemics.

An analysis of effects of heterozygosity in dairy cattle for bovine tuberculosis resistance.

Genetic selection of cattle more resistant to bovine tuberculosis (bTB) may offer a complementary control strategy. Hypothesising underlying non-additive genetic variation, we present an approach using genome-wide high density markers to identify genomic loci with dominance effects on bTB resistance and to test previously published regions with heterozygote advantage in bTB. Our data comprised 1151 Holstein-Friesian cows from Northern Ireland, confirmed bTB cases and controls, genotyped with the 700K Illumina BeadChip. Genome-wide markers were tested for associations between heterozygosity and bTB status using marker-based relationships. Results were tested for robustness against genetic structure, and the genotypic frequencies of a significant locus were tested for departures from Hardy-Weinberg equilibrium. Genomic regions identified in our study and in previous publications were tested for dominance effects. Genotypic effects were estimated through ASReml mixed models. A SNP (rs43032684) on chromosome 6 was significant at the chromosome-wide level, explaining 1.7% of the phenotypic variance. In the controls, there were fewer heterozygotes for rs43032684 (P < 0.01) with the genotypic values suggesting that heterozygosity confers a heterozygote disadvantage. The region surrounding rs43032684 had a significant dominance effect (P < 0.01). SNP rs43032684 resides within a pseudogene with a parental gene involved in macrophage response to infection and within a copy-number-variation region previously associated with nematode resistance. No dominance effect was found for the region on chromosome 11, as indicated by a previous candidate region bTB study. These findings require further validation with large-scale data.

Maternal, dominance and additive genetic effects in Nile tilapia; influence on growth, fillet yield and body size traits.

There are only few studies of dominance effects in non-inbred aquaculture species, since commonly used mating designs often have low power to separate dominance, maternal and common environmental effects. Here, a factorial design with reciprocal cross, common rearing of eggs and subsequent lifecycle stages and pedigree assignment using DNA microsatellites was used to separate these effects and estimate dominance (d2) and maternal (m2) ratios in Nile tilapia for six commercial traits. The study included observations on 2524 offspring from 155 full-sib families. Substantial contributions of dominance were observed (P < 0.05) for body depth (BD) and body weight at harvest (BWH) with estimates of d2 = 0.27 (s.e. 0.09) and 0.23 (s.e. 0.09), respectively in the current breeding population. In addition the study found maternal variance (P < 0.05) for BD, BWH, body thickness and fillet weight explaining ~10% of the observed phenotypic variance. For fillet yield (FY) and body length (BL), no evidence was found for either maternal or dominance variance. For traits exhibiting maternal variance, including this effect in evaluations caused substantial re-ranking of selection candidates, but the impact of including dominance effects was notably less. Breeding schemes may benefit from utilising maternal variance in increasing accuracy of evaluations, reducing bias, and developing new lines, but the utilisation of the dominance variance may require further refinement of parameter estimates.

Fine-mapping host genetic variation underlying outcomes to Mycobacterium bovis infection in dairy cows.

BACKGROUND: Susceptibility to Mycobacterium bovis infection in cattle is governed in part by host genetics. However, cattle diagnosed as infected with M. bovis display varying signs of pathology. The variation in host response to infection could represent a continuum since time of exposure or distinct outcomes due to differing pathogen handling. The relationships between host genetics and variation in host response and pathological sequelae following M. bovis infection were explored by genotyping 1966 Holstein-Friesian dairy cows at 538,231 SNPs with three distinct phenotypes. These were: single intradermal cervical comparative tuberculin (SICCT) test positives with visible lesions (VLs), SICCT-positives with undetected visible lesions (NVLs) and matched controls SICCT-negative on multiple occasions. RESULTS: Regional heritability mapping identified three loci associated with the NVL phenotype on chromosomes 17, 22 and 23, distinct to the region on chromosome 13 associated with the VL phenotype. The region on chromosome 23 was at genome-wide significance and candidate genes overlapping the mapped window included members of the bovine leukocyte antigen class IIb region, a complex known for its role in immunity and disease resistance. Chromosome heritability analysis attributed variance to six and thirteen chromosomes for the VL and NVL phenotypes, respectively, and four of these chromosomes were found to explain a proportion of the phenotypic variation for both the VL and NVL phenotype. By grouping the M. bovis outcomes (VLs and NVLs) variance was attributed to nine chromosomes. When contrasting the two M. bovis infection outcomes (VLs vs NVLs) nine chromosomes were found to harbour heritable variation. Regardless of the case phenotype under investigation, chromosome heritability did not exceed 8% indicating that the genetic control of bTB resistance consists of variants of small to moderate effect situated across many chromosomes of the bovine genome. CONCLUSIONS: These findings suggest the host genetics of M. bovis infection outcomes is governed by distinct and overlapping genetic variants. Thus, variation in the pathology of M. bovis infected cattle may be partly genetically determined and indicative of different host responses or pathogen handling. There may be at least three distinct outcomes following M. bovis exposure in dairy cattle: resistance to infection, infection resulting in pathology or no detectable pathology.

Cow genotyping strategies for genomic selection in a small dairy cattle population.

This study compares how different cow genotyping strategies increase the accuracy of genomic estimated breeding values (EBV) in dairy cattle breeds with low numbers. In these breeds, few sires have progeny records, and genotyping cows can improve the accuracy of genomic EBV. The Guernsey breed is a small dairy cattle breed with approximately 14,000 recorded individuals worldwide. Predictions of phenotypes of milk yield, fat yield, protein yield, and calving interval were made for Guernsey cows from England and Guernsey Island using genomic EBV, with training sets including 197 de-regressed proofs of genotyped bulls, with cows selected from among 1,440 genotyped cows using different genotyping strategies. Accuracies of predictions were tested using 10-fold cross-validation among the cows. Genomic EBV were predicted using 4 different methods: (1) pedigree BLUP, (2) genomic BLUP using only bulls, (3) univariate genomic BLUP using bulls and cows, and (4) bivariate genomic BLUP. Genotyping cows with phenotypes and using their data for the prediction of single nucleotide polymorphism effects increased the correlation between genomic EBV and phenotypes compared with using only bulls by 0.163±0.022 for milk yield, 0.111±0.021 for fat yield, and 0.113±0.018 for protein yield; a decrease of 0.014±0.010 for calving interval from a low base was the only exception. Genetic correlation between phenotypes from bulls and cows were approximately 0.6 for all yield traits and significantly different from 1. Only a very small change occurred in correlation between genomic EBV and phenotypes when using the bivariate model. It was always better to genotype all the cows, but when only half of the cows were genotyped, a divergent selection strategy was better compared with the random or directional selection approach. Divergent selection of 30% of the cows remained superior for the yield traits in 8 of 10 folds.

Genetic evaluation for bovine tuberculosis resistance in dairy cattle.

Genetic evaluations for resistance to bovine tuberculosis (bTB) were calculated based on British national data including individual animal tuberculin skin test results, postmortem examination (presence of bTB lesions and bacteriological culture for Mycobacterium bovis), animal movement and location information, production history, and pedigree records. Holstein cows with identified sires in herds with bTB breakdowns (new herd incidents) occurring between the years 2000 and 2014 were considered. In the first instance, cows with a positive reaction to the skin test and a positive postmortem examination were defined as infected. Values of 0 and 1 were assigned to healthy and infected animal records, respectively. Data were analyzed with mixed models. Linear and logit function heritability estimates were 0.092 and 0.172, respectively. In subsequent analyses, breakdowns were split into 2-mo intervals to better model time of exposure and infection in the contemporary group. Intervals with at least one infected individual were retained and multiple intervals within the same breakdown were included. Healthy animal records were assigned values of 0, and infected records a value of 1 in the interval of infection and values reflecting a diminishing probability of infection in the preceding intervals. Heritability and repeatability estimates were 0.115 and 0.699, respectively. Reliabilities and across time stability of the genetic evaluation were improved with the interval model. Subsequently, 2 more definitions of "infected" were analyzed with the interval model: (1) all positive skin test reactors regardless of postmortem examination, and (2) all positive skin test reactors plus nonreactors with positive postmortem examination. Estimated heritability was 0.085 and 0.089, respectively; corresponding repeatability estimates were 0.701 and 0.697. Genetic evaluation reliabilities and across time stability did not change. Correlations of genetic evaluations for bTB with other traits in the current breeding goal were mostly not different from zero. Correlation with the UK Profitable Lifetime Index was moderate, significant, and favorable. Results demonstrated the feasibility of a national genetic evaluation for bTB resistance. Selection for enhanced resistance will have a positive effect on profitability and no antagonistic effects on current breeding goal traits. Official genetic evaluations are now based on the interval model and the last bTB trait definition.

Genomic evaluation, breed identification, and population structure of Guernsey cattle in North America, Great Britain, and the Isle of Guernsey.

As of December 2015, 2,376 Guernsey bulls and cows had genotypes from collaboration between the United States, Canada, the United Kingdom, and the Isle of Guernsey. Of those, 439 bulls and 504 cows had traditional US evaluations, which provided sufficient data to justify investigation of the possible benefits of genomic evaluation for the Guernsey breed. Evaluation accuracy was assessed using a traditional 4-yr cutoff study. Twenty-two traits were analyzed (5 yield traits, 3 functional traits, and 14 conformation traits). Mean reliability gain over that for parent average was 16.8 percentage points across traits, which compares with 8.2, 18.5, 20.0, and 32.6 percentage points reported for Ayrshires, Brown Swiss, Jerseys, and Holsteins, respectively. Highest Guernsey reliability gains were for rump width (44.5 percentage points) and dairy form (40.5 percentage points); lowest gains were for teat length (1.9 percentage points) and rear legs (side view) (2.3 percentage points). Slight reliability losses (1.5 to 4.5 percentage points) were found for udder cleft, final score, and udder depth as well as a larger loss (13.6 percentage points) for fore udder attachment. Twenty-one single nucleotide polymorphisms were identified for Guernsey breed determination and can be used in routine genotype quality control to confirm breed and identify crossbreds. No haplotypes that affect fertility were identified from the current data set. Principal component analysis showed some divergence of US and Isle of Guernsey subpopulations. However, the overlap of US, Canadian, UK, and Isle of Guernsey subpopulations indicated the presence of gene flow, and the similarities in the subpopulations supports a common genomic evaluation system across the regions.

Efficiency of genomic prediction for boar taint reduction in Danish Landrace pigs.

Genetic selection against boar taint, which is caused by high skatole and androstenone concentrations in fat, is a more acceptable alternative than is the current practice of castration. Genomic predictors offer an opportunity to overcome the limitations of such selection caused by the phenotype being expressed only in males at slaughter, and this study evaluated different approaches to obtain such predictors. Samples from 1000 pigs were included in a design which was dominated by 421 sib pairs, each pair having one animal with high and one with low skatole concentration (≥0.3 µg/g). All samples were measured for both skatole and androstenone and genotyped using the Illumina SNP60 porcine BeadChip for 62 153 single nucleotide polymorphisms. The accuracy of predicting phenotypes was assessed by cross-validation using six different genomic evaluation methods: genomic best linear unbiased prediction (GBLUP) and five Bayesian regression methods. In addition, this was compared to the accuracy of predictions using only QTL that showed genome-wide significance. The range of accuracies obtained by different prediction methods was narrow for androstenone, between 0.29 (Bayes Lasso) and 0.31 (Bayes B), and wider for skatole, between 0.21 (GBLUP) and 0.26 (Bayes SSVS). Relative accuracies, corrected for h(2) , were 0.54-0.56 and 0.75-0.94 for androstenone and skatole respectively. The whole-genome evaluation methods gave greater accuracy than using only the QTL detected in the data. The results demonstrate that GBLUP for androstenone is the simplest genomic technology to implement and was also close to the most accurate method. More specialised models may be preferable for skatole.

Genetic contributions and their optimization.

Genetic contributions were first formalized in 1958 by James and McBride (Journal of Genetics, 56, 55-62) and have since been shown to provide a unifying framework for theories of gain and inbreeding. As such they have underpinned the development of methods that provide the most effective combination of maximizing gain whilst managing inbreeding and loss of genetic variation. It is shown how this optimum contribution technology can be developed from theory and adapted to provide practical selection protocols for a wide variety of situations including overlapping generations and multistage selection. The natural development of the theory to incorporate genomic selection and genomic control of inbreeding is also shown.

Genome-wide association study identifies novel loci associated with resistance to bovine tuberculosis.

Tuberculosis (TB) caused by Mycobacterium bovis is a re-emerging disease of livestock that is of major economic importance worldwide, as well as being a zoonotic risk. There is significant heritability for host resistance to bovine TB (bTB) in dairy cattle. To identify resistance loci for bTB, we undertook a genome-wide association study in female Holstein-Friesian cattle with 592 cases and 559 age-matched controls from case herds. Cases and controls were categorised into distinct phenotypes: skin test and lesion positive vs skin test negative on multiple occasions, respectively. These animals were genotyped with the Illumina BovineHD 700K BeadChip. Genome-wide rapid association using linear and logistic mixed models and regression (GRAMMAR), regional heritability mapping (RHM) and haplotype-sharing analysis identified two novel resistance loci that attained chromosome-wise significance, protein tyrosine phosphatase receptor T (PTPRT; P=4.8 × 10(-7)) and myosin IIIB (MYO3B; P=5.4 × 10(-6)). We estimated that 21% of the phenotypic variance in TB resistance could be explained by all of the informative single-nucleotide polymorphisms, of which the region encompassing the PTPRT gene accounted for 6.2% of the variance and a further 3.6% was associated with a putative copy number variant in MYO3B. The results from this study add to our understanding of variation in host control of infection and suggest that genetic marker-based selection for resistance to bTB has the potential to make a significant contribution to bTB control.

Assessing the impact of genomic selection against hip dysplasia in the Labrador Retriever dog.

Many purebred dogs exhibit a higher prevalence of inherited diseases compared with non-purebred dogs. One of the most popular breeds in the UK is the Labrador Retriever, which has a high prevalence of hip dysplasia resulting in high costs for surgical operations and impaired animal welfare. Considering the many complications of highly managed populations, mainly due to breeder's conventions and the resulting population structure, is of great importance for the proper development of a strategy against the disease. In this study, we have compared the utilities and performances of both genomic and phenotypic selection against hip dysplasia in a simulated population with the characteristics of the British Veterinary Association and Kennel Club (BV /KC) hip dysplasia scheme. The results confirm the potential benefits of genomic selection by showing a moderate increase of 1.15-fold (assuming a realistic accuracy of r(2) = 0.5) in response to selection due to the higher accuracy (between 0.96- and 1.32-fold, considering 0.35 ≤ r(2) ≤ 0.7) and more than a threefold increase when all the offspring in each litter are tested (between 3.25- and 4.55-fold, again considering 0.35 ≤ r(2) ≤ 0.7).

Genomic selection using beef commercial carcass phenotypes.

In this study, an industry terminal breeding goal was used in a deterministic simulation, using selection index methodology, to predict genetic gain in a beef population modelled on the UK pedigree Limousin, when using genomic selection (GS) and incorporating phenotype information from novel commercial carcass traits. The effect of genotype-environment interaction was investigated by including the model variations of the genetic correlation between purebred and commercial cross-bred performance (ρX). Three genomic scenarios were considered: (1) genomic breeding values (GBV)+estimated breeding values (EBV) for existing selection traits; (2) GBV for three novel commercial carcass traits+EBV in existing traits; and (3) GBV for novel and existing traits plus EBV for existing traits. Each of the three scenarios was simulated for a range of training population (TP) sizes and with three values of ρX. Scenarios 2 and 3 predicted substantially higher percentage increases over current selection than Scenario 1. A TP of 2000 sires, each with 20 commercial progeny with carcass phenotypes, and assuming a ρX of 0.7, is predicted to increase gain by 40% over current selection in Scenario 3. The percentage increase in gain over current selection increased with decreasing ρX; however, the effect of varying ρX was reduced at high TP sizes for Scenarios 2 and 3. A further non-genomic scenario (4) was considered simulating a conventional population-wide progeny test using EBV only. With 20 commercial cross-bred progenies per sire, similar gain was predicted to Scenario 3 with TP=5000 and ρX=1.0. The range of increases in genetic gain predicted for terminal traits when using GS are of similar magnitude to those observed after the implementation of BLUP technology in the United Kingdom. It is concluded that implementation of GS in a terminal sire breeding goal, using purebred phenotypes alone, will be sub-optimal compared with the inclusion of novel commercial carcass phenotypes in genomic evaluations.

Direct and maternal genetic relationships between calving ease, gestation length, milk production, fertility, type, and lifespan of Holstein-Friesian primiparous cows.

As the emphasis in cattle breeding is shifting from traits that increase income toward traits that reduce costs, national breeding indices are expanding to include functional traits such as calving ease (CE). However, one issue is the lack of knowledge of genetic relationships between CE and other dairy traits. The same can be said about gestation length (GL), a potential novel selection trait with considerable heritabilities and possible genetic relationships with the calving process. This study aimed to estimate the genetic relationships between CE, GL, and other dairy traits of interest using a national data set of 31,053 primiparous cow performance records, as well as to separate direct and maternal genetic effects. Chosen dairy traits included fertility (calving interval, days to first service, nonreturn rate after 56 d, number of inseminations per conception), milk production (milk yield at d 110 in milk, accumulated 305-d milk yield, accumulated 305-d fat yield, accumulated 305-d protein yield), type (udder depth, chest width, rump width, rump angle, mammary composition, stature, body depth), and lifespan traits (functional days of productive life). To allow the separation of direct and maternal genetic effects, a random sire of the calf effect was included in the multi-trait linear trivariate sire models fitted using ASReml. Significant results showed that easily born individuals were genetically prone to high milk yield and reduced fertility in first lactation. Difficult calving primiparous cows were likely associated with being high-producing, wide and deep animals, with a reduced ability to subsequently conceive. Individuals that were born relatively early were associated with good genetic merit for milk production. Finally, individuals carrying their offspring longer were genetically associated with being wide and large animals that were themselves born relatively early. The study shows that it is feasible and valuable to separate direct and maternal effects when estimating genetic correlations between calving and other dairy traits. Furthermore, gestation length is best used as an indicator trait for lowly heritable calving traits, rather than as a novel selection trait. As estimated direct and maternal genetic correlations differ, we can conclude that genetic relationships between CE, GL, and traits of interest are present, but caution is required if these traits are implemented in national breeding indices.

Estimating variance components and predicting breeding values for eventing disciplines and grades in sport horses.

Eventing competitions in Great Britain (GB) comprise three disciplines, each split into four grades, yielding 12 discipline-grade traits. As there is a demand for tools to estimate (co)variance matrices with a large number of traits, the aim of this work was to investigate different methods to produce large (co)variance matrices using GB eventing data. Data from 1999 to 2008 were used and penalty points were converted to normal scores. A sire model was utilised to estimate fixed effects of gender, age and class, and random effects of sire, horse and rider. Three methods were used to estimate (co)variance matrices. Method 1 used a method based on Gibbs sampling and data augmentation and imputation. Methods 2a and 2b combined sub-matrices from bivariate analyses; one took samples from a multivariate Normal distribution defined by the covariance matrix from each bivariate analysis, then analysed these data in a 12-trait multivariate analysis; the other replaced negative eigenvalues in the matrix with positive values to obtain a positive definite (co)variance matrix. A formal comparison of models could not be conducted; however, estimates from all methods, particularly Methods 2a/2b, were in reasonable agreement. The computational requirements of Method 1 were much less compared with Methods 2a or 2b. Method 2a heritability estimates were as follows: for dressage 7.2% to 9.0%, for show jumping 8.9% to 16.2% and for cross-country 1.3% to 1.4%. Method 1 heritability estimates were higher for the advanced grades, particularly for dressage (17.1%) and show jumping (22.6%). Irrespective of the model, genetic correlations between grades, for dressage and show jumping, were positive, high and significant, ranging from 0.59 to 0.99 for Method 2a and 0.78 to 0.95 for Method 1. For cross-country, using Method 2a, genetic correlations were only significant between novice and pre-novice (0.75); however, using Method 1 estimates were all significant and low to moderate (0.36 to 0.70). Between-discipline correlations were all low and of mixed sign. All methods produced positive definite 12 × 12 (co)variance matrices, suitable for the prediction of breeding values. Method 1 benefits from much reduced computational requirements, and by performing a true multivariate analysis.

Bias, accuracy, and impact of indirect genetic effects in infectious diseases.

Selection for improved host response to infectious disease offers a desirable alternative to chemical treatment but has proven difficult in practice, due to low heritability estimates of disease traits. Disease data from field studies is often binary, indicating whether an individual has become infected or not following exposure to an infectious disease. Numerous studies have shown that from this data one can infer genetic variation in individuals' underlying susceptibility. In a previous study, we showed that with an indirect genetic effect (IGE) model it is possible to capture some genetic variation in infectivity, if present, as well as in susceptibility. Infectivity is the propensity of transmitting infection upon contact with a susceptible individual. It is an important factor determining the severity of an epidemic. However, there are severe shortcomings with the Standard IGE models as they do not accommodate the dynamic nature of disease data. Here we adjust the Standard IGE model to (1) make expression of infectivity dependent on the individuals' disease status (Case Model) and (2) to include timing of infection (Case-ordered Model). The models are evaluated by comparing impact of selection, bias, and accuracy of each model using simulated binary disease data. These were generated for populations with known variation in susceptibility and infectivity thus allowing comparisons between estimated and true breeding values. Overall the Case Model provided better estimates for host genetic susceptibility and infectivity compared to the Standard Model in terms of bias, impact, and accuracy. Furthermore, these estimates were strongly influenced by epidemiological characteristics. However, surprisingly, the Case-Ordered model performed considerably worse than the Standard and the Case Models, pointing toward limitations in incorporating disease dynamics into conventional variance component estimation methodology and software used in animal breeding.

Estimation of historical effective population size using linkage disequilibria with marker data.

Theory hypothesizes that the rate of decline in linkage disequilibrium (LD) as a function of distance between markers, measured by r(2), can be used to estimate effective population size (N(e)) and how it varies over time. The development of high-density genotyping makes feasible the application of this theory and has provided an impetus to improve predictions. This study considers the impact of several developments on the estimation of N(e) using both simulated and equine high-density single-nucleotide polymorphism data, when N(e) is assumed to be constant a priori and when it is not. In all models, estimates of N(e) were highly sensitive to thresholds imposed upon minor allele frequency (MAF) and to a priori assumptions on the expected r(2) for adjacent markers. Where constant N(e) was assumed a priori, then estimates with the lowest mean square error were obtained with MAF thresholds between 0.05 and 0.10, adjustment of r(2) for finite sample size, estimation of a [the limit for r(2) as recombination frequency (c) approaches 0] and relating N(e) to c (1 - c/2). The findings for predicting N(e) from models allowing variable N(e) were much less clear, apart from the desirability of correcting for finite sample size, and the lack of consistency in estimating recent N(e) (<7 generations) where estimates use data with large c. The theoretical conflicts over how estimation should proceed and uncertainty over where predictions might be expected to fit well suggest that the estimation of N(e) when it varies be carried out with extreme caution.

The identification of SNPs with indeterminate positions using the Equine SNP50 BeadChip.

We have used linkage disequilibrium (LD) to identify single nucleotide polymorphisms (SNPs) on the Illumina Equine SNP50 BeadChip, which may be incorrectly positioned on the genome map. A total of 1201 Thoroughbred horses were genotyped using the Illumina Equine SNP50 BeadChip. LD was evaluated in a pairwise fashion between all autosomal SNPs, both within and across chromosomes. Filters were then applied to the data, firstly to identify SNPs that may have been mapped to the wrong chromosome and secondly to identify SNPs that may have been incorrectly positioned within chromosomes. We identified a single SNP on ECA28, which showed low LD with neighbouring SNPs but considerable LD with a group of SNPs on ECA10. Furthermore, a cluster of SNPs on ECA5 showed unusually low LD with surrounding SNPs. A total of 39 SNPs met the criteria for unusual within-chromosome LD. The results of this study indicate that some SNPs may be misplaced. This finding is significant, as misplaced SNPs may lead to difficulties in the application of genomic methods, such as homozygosity mapping, for which SNP order is important.

The unified approach to the use of genomic and pedigree information in genomic evaluations revisited.

Previous proposals for a unified approach for amalgamating information from animals with or without genotypes have combined the numerator relationship matrix A with the genomic relationship G estimated from the markers. These approaches have resulted in biased genomic EBV (GEBV), and methodology was developed to overcome these problems. Firstly, a relationship matrix, G(FG) , based on linkage analysis was derived using the same base population as A, which (i) utilizes the genomic information on the same scale as the pedigree information and (ii) permits the regression coefficients used to propagate the genomic data from the genotyped to ungenotyped individuals to be calculated in the light of the genomic information, rather than ignoring it. Secondly, the elements of G were regressed back towards their expected values in the A matrix to allow for their estimation errors. These developments were combined in a methodology LDLAb and tested on simulated populations where either parents were phenotyped and offspring genotyped or vice versa. The LDLAb method was demonstrated to be a unified approach that maximized accuracy of GEBV compared to previous methodologies and removed the bias in the GEBV. Although LDLAb is computationally much more demanding than MLAC, it demonstrates how to make best use the marker information and also shows the computational problems that need to be solved in the future to make best use of the marker data.

Phenotypic effects of calving ease on the subsequent fertility and milk production of dam and calf in UK Holstein-Friesian heifers.

The effect of calving ease on the fertility and production performance of both dam and calf was studied in approximately 50,000 and 10,000 UK Holstein-Friesian heifers and heifer calves, respectively. The first objective of this study was to estimate the effect of a difficult calving on the subsequent first-lactation milk production by estimating lactation curves using cubic splines. This methodology allows the estimation of daily milk, protein, and fat yields following calvings of differing degrees of difficulty. Losses in milk yield after a difficult calving have been quantified previously; however, estimates are generally restricted to the accumulated yields at specific days in lactation. By fitting cubic splines, gaps (in which the shape of the lactation curve can be merely guessed) between estimations were avoided. The second objective of this study was to estimate the effect of a difficult birth on the subsequent performance of the calf as an adult animal. Even though the calving process is known to involve cooperation between dam and calf, the effect of a difficult calving has, until now, only been estimated for the subsequent performance of the dam. Addressing the effects of a difficult birth on the adult calf strengthens the importance of calving ease as a selection trait because it suggests that the benefit of genetic improvement may currently be underestimated. The effect of calving ease on the subsequent reproductive performance of dam and calf was analyzed using linear regression and with calving ease score fitted as a fixed effect. Dams with veterinary-assisted calvings required 0.7 more services to conception and 8 more days to first service and experienced a 28-d longer calving interval in first lactation compared with dams that were not assisted at calving. Effects of calving ease on the reproductive performance of the adult calf in first lactation were not detected. Losses in milk yield of the dam were significant between d 9 to 90 in milk subsequent to a veterinary-assisted calving, creating a loss of approximately 2 kg of milk per day, compared with a nonassisted calving. Calves being born with difficulties showed a significant reduction in milk yield in first lactation, demonstrating the lifelong effect of a difficult birth. Compared with nonassisted calves, veterinary-assisted calves showed a loss of 710 kg in accumulated 305-d milk yield, which was significant from 129 to 261 d in milk. This suggests that from birth to production, physiological effects of a bad calving are not negated. Results furthermore suggest a beneficial effect of farmer assistance at calving on the milk yield of both dam and calf, when moderate difficulties occurred.